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ABSTRACT 



' _i ' ■ Context. The dust component of the interstellar medium (ISM) has been extensively 

^ |. studied in the past decades. Late-type stars have been assumed as the main source of 

^ , dust to the ISM, but recent observations show that supernova remnants may play a 

, role on the ISM dust feedback. 

^ Aims. In this work, we study the importance of low and high mass stars, as well as 

their evolutionary phase, on the ISM dust feedback process. We also determine the 
changes on the obtained results considering different mass distribution functions and 
star formation history. 



■ Methods. We describe a semi-empirical calculation of the relative importance of each 

. star at each evolutionary phase in the dust ejection to the ISM. We compare the 

, obtained results for two stellar mass distribution functions, the classic Salpeter initial 

mass function and the present day mass function. We used the evolutionary track 
' models for each stellar mass, and the empirical mass-loss rates and dust-to-gas ratio. 

Results. We show that the relative contribution of each stellar mass depends on the 
! ' used distribution. Ejecta from massive stars represent the most important objects for 

t I ' the ISM dust replenishment using the Salpeter IMF. On the other hand, for the present 



day mass function low and intermediate mass stars are dominant. 
Conclusions. We confirm that late-type giant and supergiant stars dominate the ISM 
dust feedback in our actual Galaxy, but this may not the case of galaxies experiencing 
high star formation rates, or at high redshifts. In those cases, SNe are dominant in the 
dust feedback process. 
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1. Introduction 

RGB and AGB stars are known as the major continuous dust producers in the Universe, 
but also SN remnant s have shown fast grain growth and dusty shells of ^ 10~^ — 10^ Mq 



( Nozawa et al 



20031) . From the classical nucleation theory, the timescale for grain growth 
in the ISM can be estimated by Tg — 4:sa{fnimiVi)~^ , where a is the dust mean size, s is 
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the material density, / the sticking probabihty, the gas phase density, rrii the atom mass 
and Vi the velocity of the i-th atom to be added onto the grain surface. Considering typical 
ISM parameters, a 0.1 — 1 /im, s ^ 2 g cm 



-3 



1 cm , r ^ 10 K, and assuming an 



efficiency / — 0.1 — 1, we find r„ < 10^ yr. 



Dust par ticles are likely to be destroyed by shock waves (jPraine &: Salpeter 



McKee 



1978; 



19891 ). Grain-grain or io n-grain collis i ons w i ll lead to the shattering process, reduc 



ing or destroying dust particles. 



Jones et al 



(|1994 119961) described the shattering process 
of ISM dust induced by SN blasts, obtaining destruction timescales of ^ 4 — 6 x 10^ yr. 
Actually, as mentioned in these works, the accurate derivation of destruction timescales 
depends on the velocity of the shock waves, the frequency of SNe and the physical proper- 
ties of the ISM, as density and temperature. Hence, the destruction timescales are shorter 
than the dust growth scale, the observed dust could not b e explained by nucleation in loco, 
but had to be recently (< 10^ yr) injected into the ISM (ITielenslliggsl ). 

It is commonly suggested that the ISM dust should be mostly originated from evolved 
low and intermediate mass stars but recent observations showed t he presence of large 



quantities of dust (Md = 10^ Mr) in the early Universe (z > 5) ( 


Hughes et al. 


1998 


Archibald et allboOl 


Dunne et al. 


2003 


Bertoldi et ahlboOS 


Maiolino et al.ll2004 


). At that 



age low and intermediate mass stars were not evolved yet and, therefo r e, SNe are recognized 



as responsible for such material ejections ( 


Todini & Ferrara 


2001 


2003; 


Sugerman et al. 


2006 


Dwek et al. 


2007 


). 



Morgan fc Edmund; 



In the post-shock phase of the S N remnant, the gas is generally cool and dense enough 



to allow dust formation and growth (jFalceta-Goncalves et al 



2003 . 



20051 ). Observationally, 



it is confirmed for SN1987A, which shows a 10 ^ Mq dust shell with a dust to gas ratio 
for heavier elements ~ 0.3. The same result was obtain ed for several other galactic and 



extragalactic SNe ([Barlow et al 



2005 



Gomez et al. 



2007). Therefor e, SNe seem to play an 



important part on the ISM dust replenishment process (jDweklll998r ). However, if the short 
dust destruction timescale is taken into account, SNe could only be the main source of dust 
of the Galaxy if it presented a high star formation rates in its recent history. 

In this work we present a semi-empirical model, in which we study the role of different 
stellar mass functions in the output of dust ejected to the ISM. The model is described in 
Sect. 2, in Sect. 3 we show the main results and present a brief discussion, followed by the 
conclusions. 



2. The Model 



To obtain the total amount of dust ejected by a certain star during all its evolution we have 
to integrate the dust mass loss rates of each evolutionary phase over the evolutionary time. 
The difficulties on performing such calculation lie on the determination of the dust mass 
loss rates at each evolutionary phase of the star, as well a s its du ration. To accomplish this 
we used the evolutionary tracks given by 



Schaller et al 



(|l992l ). which take into account 



the mass loss during the stellar evolution. These numerical calculations also provided the 
abundance of heavier elements on the stellar surface at each epoch of the stellar evolution. 
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Fig. 1. Evolutionary tracks for stars with masses ranging from 1 to 120 Mq for metallicity 
of 0.02. The box shows the values of log{M). 



We used these numerical results to determine the amount of material able to be added into 
the dust grains that is being ejected from the stars. 

The total dust mass ejected by each star over its lifetime can be determined by inte- 
grating the mass loss rate at each stellar evolutionary phase, determined by: 

r-M2 /.t2 



^{m)Md{m) dt dm, 



(1) 



where $(m) is the stellar population mass distribution, M(i{m) is dust mass loss rate of a 
given star of mass m at the time t, Mi and M2 are the limits of the mass range and ti and 
t2 are the limits of the time interval. To compute the dust mass loss rates we used: 



Md(t,m) = M(t,m)/<i(t,m)x(t,m), 



(2) 



where AI is the total mass loss rate, which is obtained empirically as described in the 
following subsection, fd is the dust-to-gas fraction of elements heavier than He and x is 
the wind metallicity, of an star with mass m at an evolutio nary time t. 



In Fig. 1 we illustrate the evolutionary track given by iSchaller et al.l (|l992r ) for stars 
with masses ranging from 0.8 to 120 Mq used in this work, for an initial metallicity of 0.02, 
and the empirical stellar mass loss rates. 

Different models of homogeneous and heterogeneous condensation theories have been 
proposed but the actual dust-to-gas fraction of stellar winds is still not completely under- 
stood. The gas acceleration is responsible for the gas rarefaction at the base of the stellar 
atmosphere and, depending on the model used, may result in very different pressure and 
temperature profiles, which are critical to the derivation of the condensation rates. As an 
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approximation, we introduce a single dust-to-gas fraction for each stellar mass, depending 
on the chemical composition. The yields for du st-to-metal fraction (fd) for the massive 
stars were the same used by 



Dwek et al 



mass stars we used the data from 



( 2007) (table 2), a nd for low and intermediate 



Morgan fc Edmund; 



(|2003h (table 1). 



2.1. The stellar mass-loss rates 

The mass loss rates, for each stellar mass at each evolutionary phase, were introduced 
empirically. For low and intermediate mass stars durin g the main sequence, giants and 



de Jager et al 



(|l988f ). The wind 



supergiants evolutionary phases, we used the data from 
velocities for giant and supergiant stars are, in general, v < 200 km s~^ and dust is not 
destroyed at heliopause. At the end of their lives, intermediate m ass-stars present very high 
mass loss rates at the post-AGB's phases (~ lO^** Mq yr~^) (jWachter et al.1 120021 ) . and 
are known as progenitors of the planetary nebulae. 

Massive stars are also know to present dust in their winds, mainly at the later stages 
of their evolution. WR and LBV stars present high mass loss rates in unstable and 
clumpy winds that f avor d ust nucleation and survival. For these objects we used the data 
given by Crowther ( 19971) . Also, about 30% of these objects are in binary systems, and 



Marchenko et al 



( 2002f) found that binary systems present high dust production rates due 
to the wind-wind shock, with 4 — 6% of the wind mass condensed into grains at WR-I-OB 
systems. On the other limit, low mass protostars were also considered since their jets and 
disks are dust growth site s. Typically, these objects present a mass-loss rate of 10~* Mq 
yr-i during 10^ - 10^ yr (jMundt et al.lll987n . 



2.2. The stellar mass distributions 

To compute the total contribution for the ejected material by all stars we may define a 
mass distribution function. To simulate the distr ibution of a yo ung stellar population we 



used the typical Salpeter single component IMF (jSalpeter 

(log M) = AM-"" 
for M* > IM0, where A = 4.43 x 10"^ and x = 1.3, and: 
4- (log M) = B exp [- (log M - log M^f /2a' 



19551) 



(3) 



(4) 



for M* < IM0, where B = 0.158, = 0.079 and a = 0.69. For the old populations, 
which exhibit a significantly lower number of massive stars, we used the Present Day Mass 
Function (PDMF). For the PDMF Eqs. ^ and (g]) are still vahd but, A = 4.4 x 10"^ and 
X = 4.37 for < logM < 0.54, A = 1.5 x 10"^ and x = 3.53 for 0.54 < logM < 1.26 and 



A = 2.5x 10-^ and x = 2.11 for log A/ > 1.26 |Chabrier 



2003 : 



Elmegreen 



20041) . 



3. Results and Discussion 

If one assumes that the dust is rapidly destroyed in the ISM (r^ ^ lO^yr), it is possible to 
naively separate in two the possible scenarios of the galactic history: i- old feedback process, 
which occurred in the presence of a large number of massive stars and, ii- recent feedback 
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Fig. 2. Dust mass, relative to total (in percent), ejected to the ISM by each stellar mass 
range during its complete evolution. The solid line represents the model using a typical 
Salpeter IMF, and the dotted line represents the model for the present day mass function 
(PDMF). 



process, with a relatively larger number of low and intermediate mass stars. Therefore, 
we must study the dependence of different stellar mass distributions on the dust feedback 
process. We performed the calculation of Eqs. ([1]) and ^ computing the ejected material 
during the evolution of all stars using the IMF and PDMF distributions. 

In Fig. 2, we show the relative amount of dust ejected by each stellar mass range to 
the ISM during all its evolution, for the IMF and the PDMF. It is noticeable that the high 
mass component is the major source of dust for the Salpeter IMF distribution. These stars 
evolves rapidly and eject large amounts of heavy atoms during the SN phase, and their 
remnants evolve to form shells of few solar masses of dust particles. On the other hand, 
the result is the opposite if use the present stellar distribution of the Galaxy. There are few 
high mass stars, and as a consequence SNe are not so frequent. Using the PDMF, the low 
and intermediate mass stars, as they leave the main sequence, are the main contributors 
of dust to the ISM. Intermediate mass stars at supergiant phase have been assumed to be 
the most important source of dust but, surprisingly, the results for the PDMF show that 
stars with < 3 Mq are dominant. It can be explained by the overabundance of low 
mass stars in this distribution. 

In the present stage of the Galaxy as the dust is destroyed, or severely changed, in short 
timescales (t^ ~ 10^ yr ), the PDMF would give more accurate abundances of dust in the 
ISM. Since the dust generated from high mass stars in the past was recycled, and with the 
absence of these objects in the current population, we may conclude that low mass stars 
are the main source of the dust observed in the present stage of the Galaxy. 

In Table 1 we show the quantitative results of the calculations described above. For a 
typical IMF the high mass stars (M* > 8 Mq) are responsible for 68% of the dust that 
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Table 1. Relative and absolute dust ISM feed- 
back contribution for different stellar component 
mass functions. 



Stellar mass range 


IMF 


PDMF 


Md (M > 8 Mq) 


68% 


1% 


Md (M < 8 Mq) 


32% 


99% 


Ma_totall{MQ) 


7 X lO'^ 


2 X 10^ 



assuming a total stellar population of 1 X 10^^ Mq. 



returns to the ISM. On the other hand, considering a PDMF, this contribution falls to less 
than 1%. To determine the absolute feedback mass, we used a total stellar population of 

1 X 10^° Mq, which gives a total dust mass Md ^ 7 x 10'' for the IMF distribution, and 

2 X 10^ for the PDMF. 

Interestingly, the results indicate that in recent starburst regions, one should expect a 
larger dust-to-gas ratio when compared to an evolved population. It could possibly be the 
reason for irregular and spiral galaxies, which present high star form ation activity, have 



more dust than the evolved elliptical galaxies (jSeaquist et al.l 12004) . However, since we 
have very different scenarios here, with different feedback and destruction timescales for 
each type of galaxy, this statement needs more detailed calculations to be tested. 

3.1. Time evolution 

In the previous calculations, to obtain the total amount of dust ejected to the ISM wc had 
to assume the mass function of the stellar population, which translates the evolutionary 
phase of the stellar component at the epoch we are studying as the dust is recycled in short 
timescales. However, in order to obtain the time evolution of the ISM dust component, we 
have to include in Eq. ([T]) the star formation rate function over history. 

We performed the evolutionary calculati o n assu ming a constant star formation rate of 



5 M0yr ^, as used bv iMorgan fc Edmunds! (|2003| ). However, differently of that work we 



took into account the dust destruction and studied its role on the total dust ejected to 
the ISM. At each time step, we calculate the total ejection from each stellar mass range 
of the current population, add new stars using the Salpeter IMF, remove stars that have 
already evolved and, finally remove the destroyed dust from the total solid component. 
In the present calculations, in order to simulate the chemical enrichment of the stars, a 
metallicity of z = O.OIzq was arbitrarily used for the initial population {t < lO^'yr), and 
z = zq for t > lO^yr. In Fig. 3 we show the results for the absolute dust mass ejected to 
the ISM. 

In both cases we can identify the dominance of massive stars on the total ejections at the 
earlier stages of the galaxy evolution. On both models, it is also noticeable the appearance of 
substantial dust amount from low and intermediate mass stars only after 10* yr. These 
objects become the major dust producers on the last 0.5 billion years on both models. 
The main differences appear on the later stages, where for considering dust destruction, 
massive stars would be responsible for 3 times less dust mass than obtained by previous 




Fig. 3. Absolute dust mass ejected to the ISM. The modeled ejection by high mass stars 
(solid thin), low and intermediate mass stars (dashed) and total (solid thick), without dust 
destruction (left panel) and considering dust destruction (right panel). 



calculations, while this proportion for low and intermediate mass stars is negligible. After 
10 billion years of evolution, we obtained a current dust mass of ^ 6 x 10^ Mq without 
dust destruction, and ~ 10'' Mq taking into account the dust destruction. 



4. Conclusions 

It is still unclear what is the main source of the dust feedback in our Galaxy. The dust must 
be formed in stars and ejected to the ISM and some authors have argued favoring cool late- 
type stellar winds, which present high mass-loss rates and dust is proved to be formed in 
these sites by observations. On the other hand, other plausible sources are the SNe ejecta. 
During the final evolutionary phases, high mass stars explode and supersonically eject a 
very rich gas to the ISM. 

Firstly, to determine the relative amount of dust ejected to the ISM for each stellar mass 
at each evolutionary phase we calculated the dust ejection during each evolutionary phase 
of the stars for different stellar mass distributions. As main result we showed that SNe are 
the main source of ISM dust feedback if a classic Salpeter IMF distribution is assumed. On 
the other hand, if we use the present day mass function, we show that the main sources 
of dust to the present Galaxy are the low and intermediate mass stars, representing more 
than 90% of the total dust mass. 

Secondly, we studied the du st feedback process along the galactic time evolution, as done 



bv [Morgan fc Edmunds! (j2003D . but including the effects of dust destruction by SN blasts. 
For simplicity we used a constant destruction rate, consistent with the current galactic 
physical parameters. During previous ages the SNe ejecta are dominant, in agreement with 
previous works. We showed that, considering the dust destruction, low and intermediate 
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mass stars are dominant for a galactic age oft > lO^yr, in a much higher proportion. The 
total dust mass of ~ lO'' M© is obtained for a star formation rate of 5 MQyr"^. 

The dust destruction timescale depends on the SNe frequency, as well as the ISM 

density and temperature. It is probable that the destruction rate was higher earlier during 
the galaetie evolution. In this case, the presented conclusions will stand, and the role of 
low and intermediate mass stars in later stages will be even higher. 
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